AC-driven quantum spins: resonant enhancement of transverse DC magnetization 
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We consider s — 1/2 spins in the presence of a constant magnetic field in ^-direction and an AC 
magnetic field in the x — z plane. A nonzero DC magnetization component in y direction is a result 
of broken symmetries. A pairwise interaction between two spins is shown to resonantly increase the 
induced magnetization by one order of magnitude. We discuss the mechanism of this enhancement, 
which is due to additional avoided crossings in the level structure of the system. 
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I. INTRODUCTION 

Magnetic resonance effects have been iQ-ngi studied 
and used for various spectroscopic methodsEftl Time- 
dependent external fields being an integral part of the 
physical setups, become also very interesting in their ac- 
tion when entering the field of spintronics and processing 
of quantum information^. While the linear response of a 
system (probe) to external fields may be typically under- 
stood in a rather complete way, nonlinear and perhaps 
nonadiabatic response effects are usually much harder to 
be analyzed. Some studies have been devoted to the ef- 
fect of second-harmonic generation. Much less is known 
about zeroth-harmonic generation, i.e. about the induc- 
tion of static fields by pure ac fields. The understanding 
of such processes is very important in any application 
involving ac fields, as additional induced dc fields may 
strongly influence the desired effects. 

The principal possibility of the zeroth-harmonic gener- 
ation can be studied using symmetry considerationstj,tl 
Analytical approaches which may also estimate the mag- 
nitude of the induced fields are typically restricted to 
various perturbation techniques with a suitable small pa- 
rameter involved. Two of us have recently applied these 
methods to the simplest cases of a spin s = 1/2 cou- 
pled to an environment (heat bath) in the presence of 
an external magnetic dc field (in z-direction)Q. While 
an additional ac field applied in ^-direction leads to the 
standard spin resonance setup used e.g. for nuclear spins, 
a tilt of the ac field direction in the x — z-plane leads to a 
broken symmetry and to a generation of a nonzero mag- 
netic moment in y-direction, i.e. perpendicular to the 
applied fields. Similar results have been also obtained 
for a spin s = 1 where the external dc field in z-directioii 
is replaced by a magnetic ion anisotropy in z-directionQ. 

This work is devoted to the analysis of spin-spin in- 
teraction and its action on the above described effects. 
The first goal is to study the generalization of the sym- 
metry considerations in the presence of an interaction. 
The second goal is to demonstrate that the interaction 
is capable of increasing the induced field strength by an 
order of magnitude, and to give an explanation for this 
result. 

The paper is organized as follows. In the next section 



we introduce the density matrix approach. Section III is 
devoted to a reconsideration of the single spin problem. 
Section IV contains the central results on two interacting 
spins. Section V concludes the paper. 



II. THE DENSITY MATRIX APPROACH 

In this section we introduce the methods we are using. 
Let us assume that a single spin or a system of interact- 
ing spins is represented by a time-dependent Hamiltonian 
H{t) — Ho + Hi(t) which includes both applied dc and 
ac magnetic fields. Here the static part is completely in- 
cluded in Hq and the action of the ac fields is described 
by a time-periodic term Hi (t) — H\ (t + T) with period 
T . In order to describe the time-dependent statistical 
evolution of such a system we use the quantum Liouville 
equation for the density matrix p(t) 



% = i[H,p] - v{p - Pe) 



(1) 



where [A, B] = AB - BA. While the specifics of the 
internal dynamics of the system are described by the first 
term on the rhs of (|l|), the second term accounts for the 
interaction of our spin system with an environment. We 
will handle this interaction in a rather qualitative way. 
The principal request this coupling should fulfill is that 
its action should preserve all symmetries of H(t) except 
time reversal symmetry. A second requirement is that 
the solution of ([|) in the absence of ac fields H±(t) = 
should be the Gibbs distribution 



Pp(H ) 



1 

2o' 



-0 H 



, Z Q = Tr(e^ H °) 



(2) 



Finally we request that all solutions of (|l|) should ap- 
proach a unique time-dependent state. 

The choice p e = pp(Ho) satisfies the above criteria. 
As shown inO for s = 1/2 it leads to equations practically 
identical with the well known Bloch equations. However 
e.g. for adiabatically slow ac fields a much better choice 
would be p e ee pp(H(t)) = ie'^W. We will discuss 
the subtle differences below. In general however we may 
state here that any of the two choices fulfills the above 
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criteria. Note that the third criterium (unique asymp- 
totic behaviour) is guaranteed, as a linear equation like 
(||) exhibits one and only one attractor, and the whole 
phase space serves as a region of attraction for this single 
attractor. Note that due to Trpp = 1 any choice with 
Trp(0) = 1 implies Trp(t) = 1 for all t. 

The value A(t) of an observable corresponding to the 
operator A is defined by 



A(t) = Tr(Ap(t)) 
and its time average is abbrevated by 
1 



A = 



T 



A(t) dt. 



(3) 



(4) 



III. REVISITING A SINGLE SPIN S = § 

We consider one spin s — 1/2 in the presence of a 
constant magnetic field in z-direction and a time-periodic 
magnetic field with angle <p to the z-axis in z — x-plane. 
The Hamiltonian is given by 



H(t) = h S z + h(t)(aS x +jS z ), 



(•5) 



where a = sin </> and 7 = cos </>. We assume that the field 
h(t) = h(t + T) has zero mean J h(t) dt = 0. 

The spin component operators are given by the Pauli 
matrices: S x>VtZ — \cr x ^y. z . The density matrix p = R+il 
can be decomposed into a real symmetric 2x2 matrix R 
and a real antisymmetric 2x2 matrix / and has three 
independent variables. Using the relation (^|), we rewrite 
the system ([!]) in terms of S x ,y,z 



S x = (h + jh(t))S y - vS x 



ah(t)S z - (h Q + jh(t))S x - vS y 
-ah(t)S y - v(S z + C) 



(G) 



where C = | tanh(/3/i /2). The structure of these equa- 
tions is equivalent to the wellknown Bloch equations pro- 
vided the relaxation times there would be chosen to be 
identicalrl Note that the symmetries of the equations (|^) 
are identical with the ones of the corresponding Bloch 
equations. 

Because of the presence of a time-periodic field, for 
times t 3> 1/v all components S XtVtZ will be also time 
periodic with period T or frequency lo — 2it/T. We can 
expand the solution of system (^|) in a Fourier series 



S Xz y, z — Aq x,y.z ~\~ 



A, 



After averaging over the period T, we have 



S X / 



Aq x; 



(7) 



(8) 



The symmetry, analysis of the equations (^) yields the 
following resultso (we considered only operations which 
conserve Trp): 

1) v = 0, 7 = , h(-t) = -h(t): 



S X * S X , Sy > Sy , S Z ► S Z , t ■ 

2) v = 0, for all 7 , h(-t) = h(t): 



Sy , S Z 



S X ^ S X ) Sy 

3) v ^ 0, 7 = 0, h(t + T/2) 



Sx 1 Si, 



-Sy , S Z 



-^S Z7 t^-i 
-h(t): 

5, , t -f t + 



T 



(9) 



(10) 



(11) 



If any of the above symmetries hold in a relevant param- 
eter case, the corresponding components of S will exactly 
vanish. E.g. for case 1) S x = 0, for case 2) S y = and for 
case 3) S x = S y = 0. The key idea now is that a proper 
choice of the parameters (including the time dependence 
of h(t)) may lead to a violation of all symmetries. As a 
consequence we expect nonzero values of the correspond- 
ing Aq X:V . This means that even if we do not apply an 
external magnetic field in y-direction, we can expect a 
nonzero S y component, which will be a function of the 
ac field frequency ix>. In order to quantify the zeroth har- 
monic generation we use the following definition of the 
strength of the effect In l ■ 



inl 



max I S y 
max I S z 



■ 100% 



(12) 



Let us emphasize here that the denominator in (|12j) is of 
nonresonant character, i.e. its value is determined by the 
maximum values of S z which is obtained for frequency 
values far from the magnetic resonance and practically 
coincides with the equilibrium value for zero ac fields. In 
contrast the enumerator gets its maximum (if nonzero et 
al) inside some resonance window in u>. 

InQ a perturbation theory was developed for strong dis- 
sipation Here we will account for the effect as- 
suming a small amplitude of the ac field h(t) = ecos(wt), 
e«l and e < v. 
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A Qz = -C 



+ 0(O 
,a 2 C(^+u; 2 + h 2 ) 
2A 



+ o(e 3 ) 



(13) 
(14) 

(15) 



where A = (hi - u 2 ) 2 + v 2 (v 2 + 2h 2 + 2uj 2 ). The de- 
pendence of Ao XiV . z on lo shows the expected resonant 
character. A nonzero component of S y , being a nonlin- 
ear response result, appears in second order of the ac field 
amplitude e. 

For the case of the magnetic resonance u> = ho we find 
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A 0: 



4>1 
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a^yC 



(16) 
(17) 



where f = 

s ho 

To test the possible magnitude of the effect, we com- 
pute the induced moments numerically following!! The 
results are the solid lines in Fig.[iJ While the resonance 
character is clearly observable, the strength Inl ~ 3.0% 
is rather small. 




FIG. 1. 
h = 3.0, 
lines: p e = 



i> = f , p = 10, v = 0.1, 
Pg(Ho). Dashed lines: p, 



Chosen parameters are: 
h(t) = V2cos(ujt). Solid 



In addition, one can see that in the adiabatic limit to — > 
S y keeps a small but nonzero value. This is in contrast 
to the expectation that for very slowly varying ac fields 
the spin system should momentarily adjust to the actual 
field value, which implies a zero average y-component 
of the induced moment. As already discussed above, a 
much better choice in this limit would be p e = pp{H{t)). 
A computation with this choice is shown also in Fig.l 
(dashed lines). We find that while the overall difference 
in the results for the two choices for p e are small, the 
adiabatic limit is correctly reproduced by the second one 
as expected. 
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FIG. 2. log S x ,y,z as function of log^ for u) = 1.5, other 
parameters as in Fig.|l]. S x : solid line, S y - dashed line, S z - 
dotted line. 

It is also instructive to study the dependence of the 
values of S x , y , z as a function of the dissipation v. In 
Fig.^ we plot this dependence for the subresonance case 
u> = 1.5 from Fig.|l|. We observe that S z shows a kink 
around v — 0.1. The S x curve provides with the a ~ 1/V 2 
decay for large v, as expected from perturbation theory. 
Except for a small dip around v = 2 it then saturates 
at a nonzero value for small v values. This is in full 
agreement with the above symmetry statements. Finally 
the S y curve shows a ~ l/i/ 3 decay for large v, again as 
expected from perturbation theory. Notably S y tends to 
zero also for small dissipation values, in full accord with 
the above symmetry statements. 



IV. TWO INTERACTING SPINS 

In order to enhance the value of Inl we may either 
consider larger spins or effects of interactions between 
spins. IrJH the case of s = 1 has been considered, with 
the result that Inl is not significantly changed as com- 
pared to s = 1/2. Consequently we turn in this section 
to the consideration of spin interactions. In the following 
we will study a system of two interacting spins s = 1/2, 
which could be realized through the interactions e.g. in 
a molecule. As we are interested in the qualitative ef- 
fects such an interaction might induce, we choose the 
mathematically simple version of exchange interaction. 
We think that e.g. dipole-dipole interactions as realized 
between nuclear spins will not significantly alter the pic- 
ture. 

Our starting Hamiltonian is given by 



H 



h (Sl + S 2 z ) 

J x S , r S „ 



JySySy 



h(t)(a(Sl+S 2 x ) 



i{sl + sl)) 



(18) 



with h(t + T) = h(t), a = sin(0) and 7 = cos(</>). For 
two spins s = 1/2 the above operators take the form 
,,2 (g) 1, Si z = 1 (g) ±a X:Z where a XtZ are Pauli 



$x,z — 
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matrices, 1 denotes the 2x2 unity matrix and (^) stands 
for the tensorial product of operators (matrices). 

In this case the matrix p = R + il contains 15 inde- 
pendent variables. 

We define the total components of spins as S XtV>z = 
k( S iyz + S l,y,z)- With the hel P of © it; follows 



-R34) 



Sy = -(fa + il3 + 124 + hi) 

S z = ~ ^ 44 ^ 



(19) 



The symmetry analysis of this system is identical to 
the case of one spin s = 1/2 (|ll|-^o|). The existence of 
additional parameters (exchange coefficients) does not re- 
duce the symmetries nor add some. The only important 
additional point is that changes of signs of the expecta- 
tion values of a given spin component operator have to be 
done simultaneously for both spins once the interaction 
is nonzero. 



A. Isotropic exchange 



Let us first consider the case of isotropic exchange J x — 

tjy J z J • 

First we note that with the help of a unitary trans- 
formation we can transform the Hamiltonian ( |l8| ) into 
a triplet-singlet representation. The singlet state (total 
spin s = 0) decouples from the three triplet states (total 
spin s = 1). The time-dependent part of (|l^ ) is contained 
in the triplet part. For a given value of the ac field the 
eigenvalues of the singlet state E and the three triplet 
states £1,2,3 are given by 



-J 



E l = --J + Va 2 + 4c 2 , 
E 3 = --J - Va 2 + 4c 2 



(20) 
(21) 



where a = ho + jcos(ojt) and c = cos(ojt) . 

On one side the singlet state is not interacting with the 
triplet states and all spin component expectation values 
in this singlet state are exactly zero. On the other side 
the singlet state is consuming statistical weight. Thus the 
dynamics of the system of two spins one-half interacting 
via isotropic exchange and coupled to a heat bath can be 
reduced to the dynamics of a single spin s = 1 (triplet 
states) with a reduced statistical weight. For the case of 
Pe = Pp(Ho) this reduced statistical weight is given by 



Tr(e- pH o) 
Tr ( e -/3i/ ) 



where H' is obtained from Hq by excluding the singlet 
state. Using Tr(e~@ H °) = 2e^~ cosh(/3/i ) + and 
Tr ( e -/3ff ) = r r ( e -/3ffo) + e 



3/3.7 



we find 



2cosh(/% ) + 1 
2cosh(/3ft, ) + 1 + e~P J 



(23) 



(22) 



The coefficient (23) tends to 1 for large ferromagnetic 
interaction J — > 00 and to for large anti-ferromagnetic 
interaction J — > — 00. 

As we know the value of Inl for J = (it is equiva- 
lent to the case of a single spin s = 1/2 discussed in the 
previous section) we may obtain the values for Inl and 
all other numbers including curves as in Fig.l using ([23]). 
First we note that antiferromagnetic interactions J < 
reduce the effect by making the singlet state energeti- 
cally favourable. This is reasonable as for zero interaction 
both spins will on average point to some direction. An- 
tiferromagnetic interaction will favour antiparallel spin 
orientations and thus reduce the induced moments. How- 
ever ferromagnetic interaction, making the singlet state 
energetically unfavourable, will enhance the effect. An 
upper bound for the maximum enhancement can be eas- 
ily obtained from (^3|). In the most favourable case the 
statistical weight of the triplet states will increase up to 
4/3 of the J — case. Thus the case of isotropic ferro- 
magnetic exchange, while leading to an enhancement of 
the induced moments, does not change the results drasti- 
cally as compared to zero interaction. We will now turn 
to anisotropic exchange interaction and see that the sit- 
uation is going to change dramatically. 



B. Anisotropic exchange 

To provide with a systematics in the search for ex- 
change values which may significantly enhance the in- 
duced magnetic moments we first mention that as for 
isotropic exchange, also the general case allows for a 
unitary transformation which separates a singlet state 
(total spin s = 0) from three triplet states (total spin 
s = 1). The singlet state will become unimportant if its 
consumed statistical weight is small. This can be typi- 
cally obtained by choosing ferromagnetic exchange inter- 
actions. Second we recall that for isotropic exchange the 
dependence of the triplet eigenenergies for a given value 
of the ac field h (frozen time, or adiabatic limit) as a 
function of this ac field value h will show up with a pair 
of states performing an avoided crossing at h — and a 
third state being independent of h. This follows immedi- 
ately from ((2^). The minimum distance between the two 
states Ei and £3 is given by the Zeeman splitting 2ho, as 
in the case of a single spin s = 1/2 (up to the factor 2). 
Also we note that the ac field term induces interactions 
and thus transitions between these two states, leading to 
the standard magnetic resonance. 

For anisotropic exchange the situation changes. The 
third level E2, which was not participating in the res- 
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onance dynamics for isotropic exchange, now gets in- 
volved in the dynamics due to additional nonzero ma- 
trix elements for anisotropic exchange. As a consequence 
we may expect resonant response at different frequen- 
cies (energies) corresponding to the different level spac- 
ings Ei — £3 and Ei — E3. Especially important is, that 
we may consider e.g. a resonant response for transitions 
from the lowest level to the highest one of the triplet 
states, while having at the same time a narrowly avoided 
crossing between the lowest state and the intermediate 
one. The idea behind this scenario is that i) small en- 
ergy differences in the vicinity of the avoided crossing 
will lead to small denominators in the equation for the 
density matrix and ii) an additional avoided crossing pro- 
vides with a strong nonlinear dependence of the wave- 
function of the lowest triplet state upon variation of h. 
So our strategy is to find exchange values which realize 
such a situation. Then we choose amplitudes of the ac 
field which will cover the avoided crossing, and temper- 
atures small enough such that even in the point of the 
avoided crossing the statistical weight of the intermediate 
state is small as compared to the triplet ground state. 

The example we will demonstrate corresponds to a fer- 
romagnetic easy plane exchange interaction J x = J y = 
5.0 and J z — 0. The dependence of the levels on the 
amplitude of the frozen ac field is shown in FigJ^. 
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frozen ac field value 
FIG. 3. Energy level dependence on frozen ac field value 
(parameters see text). Arrows indicate two possible resonant 
excitations. 
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FIG. 4. S y as function of u). Chosen parameters are: 
h = 3.0, <j> = f , p = 10.0, v = 0.1, J x = J y = 5.0 and 
J z = 0.0. In this case h(t) — v / 2cos(o»t). 



In Fig.|f| the results of numerical simulation for S y are 
shown for = 10 and v = 0.1. We do not observe a 
strong response at the expected frequency value 5.8, but 
we do observe in addition to the standard response at 
lu = ho = 3 a response at to » 0.5 which corresponds 
to the resonant excitation in the region of the narrow 
avoided crossing in Fig||. The value of Inl is already 
pretty large here, of the order of 7%. 

In order to enhance the response at w « 5.8, we de- 
crease the value for the dissipation constant by two orders 
of magnitude to v — 10~ 3 . The resulting curve is shown 
in Fig.[|. We clearly find a strong enhancement of the res- 
onance at lo w 5.8. In this case the value of Inl ~ 18%, 
which is the largest value we achieved to find. To char- 
acterize the dependence of this strong resonance on the 
parameters and v we plot different curves for additional 
values of v (Fig. ||) and (Fig. 0). The dependence on the 
inverse temperature shows the expected saturation at low 
temperatures. The dependence on v shows a surprising 
non-monotonic behaviour. While it would be tempting 
to discuss it in terms of so-called stochastic resonance, we 
think that a much simpler explanation is due to the fact 
that in the limit of zero dissipation time reversal sym- 
metry ( |Iof ) will be restored. Consequently in the limit 
v — > the values of the y-component of the induced dc 
magnetic moment will vanish for all frequencies. 



The minimum distance between the two lowest states 
in the avoided crossing region is 0.5. We now expect a 
possible strong enhancement of the induced moments in 
the frequency range of the resonance indicated by the ar- 
row u> ~ 5.8 which should be most effective if the temper- 
ature is smaller than the avoided crossing splitting, i.e. 
> 2. In addition we expect to see another resonance 
at frequencies corresponding to the smallest splitting be- 
tween the two lowest lying states, i.e. around to = 0.5. 



5 



0.08 



0.06 



.0.04 



Ivi 



0.02 



-0.02 



1 2 3 4 5 6 

CO 

FIG. 5. S y as function of u. Parameters as in Fig.[| except 
v = 10~ 3 
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FIG. 6. S y versus w in the strong resonance region w « 5.8 
for different values of v. Other parameters as in FigM 
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FIG. 7. S y versus w in the strong resonance region ui ~ 5.8 
for different values of /3. Other parameters as in Fig.[l| except 
v = 10~ 3 . 



V. CONCLUSION 

In conclusion we would like to stress once again that 
strong AC magnetic fields may generate zeroth harmonics 
in a spin system. This effect has resonant behaviour as 
a function of the AC field frequency and could be a new 
tool for the spectroscopic study of magnetic systems. 

We have shown that the induction of a dc magnetic mo- 
ment perpendicular to applied magnetic fields can be sig- 
nificantly enhanced for systems of interacting spin pairs 
as compared to single spins. This enhancement of the 
nonadiabatic nonlinear response effect is due to the pres- 
ence of additional energy levels in a spin pair system, 
which allows to obtain additional avoided crossings be- 
tween levels upon applying external fields. These avoided 
crossings strongly contribute to the nonlinear response 
through the dramatic change of wavefunctions as one 
sweeps through such a crossing. The presence of avoided 
crossings is typically associated with chaos and noninte- 
grability in corresponding classical systems. The case of 
isotropic exchange reduces the system to a single spin one 
and consequently does not lead to a strong enhancement 
of the induced moments. 

The particular case of two interacting spins could be 
realized in stable biradicals, which provides us with the 
system of two spins located on the same molecule. Ex- 
perimental possibility of observation of the above effect 
could be connected with diluted paramagnetic media at 
very low temperature, because only at low temperature 
the rate of dissipation of energy v may become very small, 
thus, preventing the heating up of the system and nar- 
rowing the width of resonance, which is very important 
for the observation. 

Other candidates for such an observation are dielectric 
glasses based on SiO x . The low temperature proper- 
ties of such glasses are connected with the existence of 
systems of double-well tunnel centers. The dynamics of 
these tunnel centers maps exactly onto the dynamics of 
spin systems. Thus, under the action of an AC electric 
field a permanent polarisability component could be in- 
duced. Again low temperatures are a key requirement 
for the observation of this effect. 
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